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Abstract 

Proceeding from optical analogy we propose a new physical interpreta- 
tion of unitarity saturation leading to antishadowing as a reflective scatter- 
ing. This interpretation of antishadowing is related to the non-perturbative 
aspects of strong interactions and follows from the specific property of the 
unitarity saturation when elastic S'-matrix S(s, b) 1 6=0 - ► — 1 at s — > 00. 
The analogy with Berry phase and experimental consequences of the pro- 
posed interpretation as reflective scattering at the LHC and in the cosmic 
rays studies are discussed. 
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Introduction 



The fundamental problems of the nonperturbtive QCD are related to confinement 
and spontaneous chiral symmetry breaking phenomena. Those phenomena deal 
with collective, coherent interactions of quarks and gluons, and results in for- 
mation of the asymptotic states, which are the colorless, experimentally observ- 
able particles. Among the problems of strong interactions, the total cross-section 
growth with energy constitutes one of the important questions. The nature of this 
energy dependence is not completely understood since underlying microscopic 
mechanism is due to the nonperturbative QCD. An essential role here belongs to 
elastic scattering where all hadron constituents interact coherently and therefore 
elastic scattering can also serve as an important tool to the confinement studies. 
Owing to the experimental efforts during recent decades it has became evident 
that the coherent processes survive at high energies. 

General principles play a guiding role in the soft hadron interaction studies, 
and, in particular, unitarity which regulates the relative strength of elastic and 
inelastic processes is the most significant one. It is important to note here that 
unitarity is formulated for the asymptotic colorless hadron on-mass shell states 
and is not directly connected to the fundamental fields of QCD — quarks and 
gluons. The same is valid for the analyticity which is relevant for the scattering 
amplitudes of the observable particles only. As it was noted in [JT]], it is not clear 
what these fundamental principles imply for the confined objects. 

Even the extension of S'-matrix and s-channel unitarity to the off -mass shell 
particle scattering leads to the significant changes in the predictions for the be- 
havior of the observables: it does not rule out now the power-like asymptotical 
behavior of the total cross-sections, i.e. unitarity in this case does not lead to the 
Froissart-Martin bound in the small £ -region 0. 

Our goal here is to provide a new physical interpretation for the scattering 
mode based on the saturation of the unitarity relation for the on-mass shell parti- 
cles, when elastic scattering not only survives, but prevails at super high energies 
in hadron collisions at small impact parameter values. We discuss in this note this 
rather unexpected behavior on the base of optical analogy, and list some of the 
respective experimental signatures in the studies of hadronic interactions at the 
LHC and in cosmic rays. It is important to note that some of these predictions are 
due to this new interpretation of unitarity saturation as a reason for the appearance 
of the reflective scattering at superhigh energies. 
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1 Saturation of unitarity and energy evolution of the 
geometric scattering picture 

Unitarity or conservation of probability, which can be written in terms of the scat- 
tering matrix as 

SS+ = 1, (1) 

implies an existence at high energies of the two scattering modes - shadowing and 
antishadowing. Existence of antishadowing has been known long ago flUHl, in 
the context of the rising total cross-sections and transition beyond the black disk 
limit it was discussed in the framework of the rational unitarization scheme in 0. 
The most important feature of antishadowing is self-damping j3]| of the inelastic 
channels contribution. Corresponding geometric picture and physical effects will 
be discussed further. As it was noted in Introduction, the operator S is defined in 
the Hilbert space spanned over vectors corresponding to the observable physical 
particles and the unitarity is formulated in terms of the physical particles. 

The attempts to construct S-matrix in perturbative QCD [9] have difficulties 
related to the problem of confinement JH [TOj and the problem of the undefined 
scale dependence entering S'-matrix through a s ifTTTl . Moreover, as it was shown 
in [fT2l . the color fields can be represented by harmonic oscillators (gluons) only 
in the limit when the strong interaction coupling constant tends to zero, i.e. at 
Q 2 -» oo. 

After these necessary remarks, one can proceed further with unitarity condi- 
tion, apply the standard procedure and derive from Eq. (OQ) unitarity relation for 
the partial wave amplitudes fi(s) : 

lmf l (s) = \f l (s)\ 2 + r ]l (s), (2) 

where elastic S'-matrix is related to the amplitude as 

S l (s) = l + 2if l {s) (3) 

and r)i(s) stands for the contribution of the intermediate inelastic channels to the 
elastic scattering with the orbital angular momentum /. The relation © turns out 
to be a quadratic equation in the case of the pure imaginary scattering amplitude. 
Therefore, the elastic amplitude appears to be not a singe- valued function of r/j. 
But, only one of the two solutions of the equation corresponding to the relation 
© is considered almost everywhere 

fi{s)= l -{l- i-e. |/j|<1/2, (4) 

while another one 

f l {s)= % -{l + s /l-±r n {s)), i.e. 1/2<|/,|<1 (5) 
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is ignored. Eq.© corresponds to the typical shadow picture when elastic ampli- 
tude fi(s) for each value of I is determined by the contribution of the inelastic 
channels r]i(s) and for small values of the function r]i(s): fi(s) ~ ii]i(s). Note, 
that the value \fi\ = 1/2 corresponds to the black disk limit where absorption 
is maximal, i.e. 7#(s) = 1/4. Unitarity limit for the partial wave amplitude is 
unity, i.e. twice as much as the black disk limit. The saturation of the unitarity 
limit is provided by Eq.©; it leads at the small values of r)i(s) to the relation 
fi(s) ~ i(l — i]i(s)). It means that elastic amplitude is increasing when contribu- 
tion of the inelastic channels is decreasing. Therefore, the term of antishadowing 
was used. 

It is important to consider the arguments in favor of the solution © neglecting. 
They are simple: it is well known that analytical properties in the complex t- 
plane lead to decrease with I at large values of / > L(sf\ of both the amplitude 
fi(s) and the contribution of inelastic channels r/z(s) at least exponentially, i.e. 
lim^oo fi(s) = and lim^oo r]i(s) = 0. It is evident that the solution © does 
not fulfill the requirement of simultaneous vanishing fi(s) and r]i(s) at I — > 00 and 
this is the reason for its neglecting. But the requirement of simultaneous vanishing 
fi(s) and rji(s) effective at large values of / only and, respectively, the solution © 
should be neglected at large values of /. At small and moderate values of I both 
solutions © and © have equal rights to exist. 

But how both the above mentioned solutions can be reconciled and realized 
in the uniform way and what physics picture can underlie the solution ©, which 
anticipates saturation of the unitarity limit for the partial wave amplitude? In 
what follows we consider the case and its consequences. To provide a geometric 
meaning to the scattering picture we will use an impact parameter representation. 
The unitarity relation written for the elastic scattering amplitude f(s,b) is similar 
to © in the high energy limit and has the following form 

lmf(s,b) = \f(s,b)\ 2 + r ] (s,b). (6) 

There is no universally accepted method to implement unitarity in high energy 
scattering |0]|. The recent approaches to solve the problem of high-energy limit 
in QCD based, e. g. on the AdS scattering, apriori suppose the eikonal function 
corresponding to the black disk [5J. The reduced unitarity limit for the amplitude 
based on Eq. © is used also in the models of gluon saturation (cf. O and 
references therein). 

In principle, a choice of particular unitarization scheme is not completely am- 
biguous. The two above mentioned solutions of unitarity equation © and © 
can be naturally reconciled and uniformly reproduced by the rational (U -matrix) 

'it should be noted that L(s) ~ ^/s In s. 
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form of unitarization. The arguments based on analytical properties of the scat- 
tering amplitude were put forward lfT3l in favor of this form. In the [/-matrix 
approach the elastic scattering matrix in the impact parameter representation is 
the following linear fractional transform: 

l+iU(s,b) 

S(s,b) = - — ( hV (V) 

1 — iU(s, b) 

U (s, b) is the generalized reaction matrix, which is considered to be an input dy- 
namical quantity. The transform © is one-to-one and easily invertible. Inelastic 
overlap function r)(s, b) can also be expressed through the function U (s, b) by the 
relation 

and the only condition to obey unitarity in the form of Eq. © is ImU (s, b) > 0. 

Another way to warrant unitarity is to represent elastic S-matrix in the expo- 
nential form: 

S(s,b) =exp[2i5{a,b)], (9) 

where 5(s,b) is a phase shift [5(s, b) = 5r(s,6) + iSi(s, b)] and the inequality 
Si(s, b) > is needed to satisfy unitarity. Both representations © and © provide 
\S(s, 6)| < 1 but contrary to ©, S ^ for any finite value of 6(s, b) when S is 
written in the form ©. To trace a further difference between them, let us consider 
for simplicity the case of pure imaginary [/-matrix and make the replacement 
U — > iU. The S'-matrix has the following form 

1 + U(s, b) 

At this point we should make an important remark on the models, which use ra- 
tional or exponential forms of the amplitude unitarization. We note that most of 
the models provide increasing dependence of these functions with energy (e.g. 
power-like one) and their exponential decrease with impact parameter b. Thus, 
any model of this kind for [/-matrix is not compatible with any similar eikonal 
model, predicting crossing the black disk limit. The value of energy correspond- 
ing to this limit for central collisions S(s, 6) | &=o = will be denoted as s and it is 
determined by the equation U (s, £») | t,=o = 1- Thus, in the energy region s < s the 
scattering in the whole range of impact parameter variation has a shadow nature 
and correspond to solution ©, the S matrix varies in the range < S(s, b) < 1. 
But when the energy is higher than this threshold value so, the scattering pic- 
ture at small values of impact parameter b < R(s) corresponds to the solution 
Eq. ©, where R(s) is the interaction radius. The S'-matrix variation region is 
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— 1 < S(s,b) < at s > s and b < R(s). We are going to discuss emerging 
physical picture of the scattering in this particular region of impact parameters 
and very high energies. The schematic energy evolution of S'-matrix of the impact 
parameter dependence is depicted on Fig. 1. It is evident that at s > sq there is 
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Figure 1: Impact parameter dependence (schematic) of the function S(s,b) for three 
values of energy si < So < S2- 

a region of impact parameters where S'-matrix is negative, i.e. phases of incom- 
ing state and outgoing state differ by n. Thus, there is a close analogy here with 
the light reflection off a dense medium, when the phase of the reflected light is 
changed by 180°. Therefore, using the optical concepts [14J, the above behavior 
of S(s, b) should be interpreted as an appearance of a reflecting ability of scatterer 
due to increase of its density beyond some critical value, corresponding to refrac- 
tion index noticeably greater than unity. In another words, the scatterer has now 
not only absorption ability (due to presence of inelastic channels), but it starts to 
be reflective at very high energies and its central part (b = 0) approaches to the 
completely reflecting limit (S = — 1) at s — > oo. It would be natural to expect that 
this reflection has a diffuse character. In another words we can describe emerging 
physical picture of very high energy scattering as scattering off the reflecting disk 
(approaching to complete reflection at the center) which is surrounded by a black 
ring. The reflection which is a result of antishadowing leads to S(s, b)\ b=0 — > — 1 
asymptotically. 

The inelastic overlap function r/(s, b) gets a peripheral impact parameter de- 
pendence in the region s > s (Fig. 2). Such a dependence is manifestation of the 
self-damping of the inelastic channels at small impact parameters. The function 
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b=R(s) 



Figure 2: Impact parameter dependence of the inelastic overlap function in the unita- 
rization scheme with antishadowing. Arrows indicate directions of the energy evolution. 



rj(s, b) reaches its maximum value at b = R(s), while the elastic scattering (due 
to reflection) occurs effectively at smaller values of the impact parameters, i.e. 
(b 2 ) el < (b 2 ) mel . Note that 

(b% = - fb 2 da^- Hb^db 2 , 
v% J °i Jo d ° 

where i = tot, el, inel and 

T n, n 1 d&tot i /■/ tm2 1 d(T e i 1 d(Ji ne i 

The quantity (b 2 ) is a measure of the particular reaction peripherality and the 
following sum rule takes place for {b 2 ) ine i. (b 2 ) inel a inel (s) = E n > 3 ( 62 }n CT n( s )' 

It is useful to return to the exponential representation for the S-matrix © at 
this point. It is evident that this form with pure imaginary phase shift (eikonal) 
would never give the negative values of the function S(s, b), it will always vary in 
the range < S(s, b) < 1. However it is not the case when 5r(s, b) is not zero. If 
5r(s, b) = 7r/2, the antishadowing can be reproduced in the exponential form of 
the unitarization, i.e. the function S(s, b) will vary in the range —1 < S(s, b) < 0. 
In order to combine shadowing at large values of b with antishadowing in central 
collisions the real part of the phase shift should have the dependence 

6 R ( s ,b) = ^6(R(s)-b). 

Such a behavior of Sr(s, b) just takes place in the [/-matrix form of unita- 
rization. Indeed, the phase shift 8(s, b) can be expressed in terms of the function 
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U(s,b) as following 

s{s , b) = W^HM (n) 

2i 1 + U(s, b) 

It is clear that in the region s > s the function 8(s, b) has a real part n/2 in the 
region < 6 < while goes to infinity at 6 = i?(s) (Fig. 3). The 
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Figure 3 : Impact parameter dependence of the real (left panel) and imaginary part ( right 
panel) of the phase shift at s > sq. 

picture leads to discontinuity in the phase shifj^ and it resembles the scattering on 
the cut-off optical potentials lfi~6l . 

It also leads to another interesting similarity, namely it allows one to consider 
8r as an analog of the geometric Berry phase in quantum mechanics which ap- 
pears as a result of a cyclic time evolution of the Hamiltonian parameters [fT71 . 
This interesting phenomenon can be observed in many physical systems [TT8l [T9l 
and is, in fact, a feature of a system that depends only on the path it evolves along 
of. In the case of pure imaginary elastic scattering amplitude the contribution of 
the inelastic channels i] can be considered as a parameter which determines due 
to unitarity (but not in a unique way) the elastic S-matrix. We can vary variable s 
(and/or b) in a way that the parameter r] (which has a peripheral 6-dependence, cf. 
Fig. 2) evolves cyclically from rji < 1/4 to r] ma x = 1/4 and again to the value r/f, 
where r/f = r/i (loop variation). As a result the non-zero phase appears (Sr = n/2 
at b < R(s)) and this phase is independent of the details of the energy evolution. 

Thus, we can summarize that the physical scattering picture beyond the black 
disk limit evolves with energy by simultaneous increase of the reflective ability 
(i.e. \S(s, b) | increases with energy and 5r = n/2), and decrease of the absorptive 
ability 1 — \S(s, b)\ 2 at the impact parameters b < R(s). 

2 The discontinuity of Si(s, b) was noticed first by V.A. Petrov lfT51 . 
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We address now the elastic scattering amplitude F(s,t), which is a Fourier- 
Bessel transform of the f(s, b). In the [/-matrix unitarization method it is deter- 
mined by the singularities in the impact parameter complex (3 = b 2 plane. Those 
singularities include the poles which positions are determined by the solutions of 
the equation 1 + U(s, /3) = and the branching point at (3 = which follows 
from the spectral representation for the function U(s,/3) 



a J to 

With account of the analytical properties dictated by Eq. (TT2l) we can use param- 
eterization in the form EOl l22l 



where M is the total mass of N constituent quarks in the colliding hadrons and £ is 
a parameter. Contribution of the poles determines the elastic scattering amplitude 
in the region of small and moderate values of —t. The amplitude dependence 
in this region provides diffraction peak and shows up dip-bump structure of the 
differential cross-section. It reproduces also Orear behavior at larger values of —t. 

The pole and cut contributions are decoupled dynamically when g(s) — > oo 
at s — > oo ll20l l22"l . At large angles the contribution from the branching point is 
a dominating one. The large angle or small impact parameter scattering being a 
result of the reflection has a power-like angular distribution dependence: 



The power-like dependence of the differential cross sections in large angle scat- 
tering closely interrelates with the rise of the total cross sections at high energies 
since both are determined by the dependence of g(s) — > oo at s — > oo in the 
unitarity saturating scheme. 

For the symmetric case of pp-interaction the scattering is the same in the for- 
ward and backward hemispheres. The more interesting case is the one where 
interacting particles are not identical, for example, 7riV-scattering. In this case 
the behavior of the differential cross-section in the forward hemisphere is com- 
pletely analogous to the above case of pp-scattering where overall behavior of the 
differential cross-section incorporates diffraction cone, Orear type and power-like 
dependencies. But in the backward hemisphere the poles contributions are sup- 
pressed compared to the cut contribution in the whole region of the variation of 
the variable u , i.e. the power-like dependence will take place at all values of u 
and there would be no diffraction cone and Orear type dependence [|20l in the 
backward hemisphere. It is not surprising, indeed, if we will recollect reflecting 
nature of the scattering at small impact parameters. 




(12) 




(13) 




(14) 
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2 Observable effects at the LHC energies 



In the [/-matrix unitarization scheme asymptotical behavior of the cross-sections 
and the ratio of elastic to total cross-section are different from the asymptotic 
equipartition of elastic and inelastic cross-sections based in the black disk limit 
when 

&el(s) / (J tot (s) -> -. 

Under the [/-matrix unitarization the rise of the elastic cross-section is predicted 
to be steeper than the rise of the inelastic cross-section beyond some threshold 
energy. Asymptotically 

crtot(s) ~ er eZ (s) ~ ln 2 s, a ine i(s) ~ Ins. (15) 

This is due to saturation of the upper unitarity limit for the partial-wave amplitude 
in the [/-matrix approach |/;| < 1, while the restriction \fi\ < 1/2 dictated by the 
black disk limit leads to the above mentioned equipartition of the cross-sections. 
Note, that despite the asymptotics for a e i and a ine i are different, the quantities 
(b 2 ) e i and (b 2 ) ine i have the same energy dependence, proportional to In 2 s at s — > 
oo. 

The above asymptotic dependencies take place for various forms of the func- 
tion U (s, b). Explicit forms^l for the function U (s, b) can be obtained, e.g. using 
geometrical, Regge or chiral quark models for the [/-matrix. Of course, it is use- 
ful to have numerical estimates for the cross-sections at the LHC energies. These 
estimates are model dependent ones and can vary rather strongly depending on the 
choice of the particular model parameterizations for [/-matrix. However, both the 
Donnachie-Landshoff and dipole Pomeron parameterizations of [/-matrix, used in 
[12711 and the model [22j are in agreement that the black disk limit will be passed 
at yfs = 2 TeV. The latter model provides the following values at the LHC energy 
y/s = 14 TeV: a to t — 230 mb and a e i/a tot — 0.67 [|23l . Thus, there is an in- 
teresting possibility that the reflective scattering mode could be discovered at the 
LHC by measuring o d / o tot ratio which would be greater than the black disk value 
1/2. However, the asymptotical regime (fT5l) is expected in the model at \fs > 100 
TeV only. 

Proceeding from an increasing weight at the LHC energies of the reflective 
scattering compared to the shadow scattering, we can suppose that the dip-bump 
structure in the da/dt in pp scattering might be less prominent at large values of 
—t and hadronic glory effect (enhancement of backscattering probability) would 
be observed. Unfortunately, it is difficult to give a more definite predictions for 

3 Of course, quantitative predictions are different for the different models while qualitative 
features coincide. 
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differential cross-sections, since other effects such as real part of the amplitude 
and/or contributions from helicity-flip amplitudes should be taken into account. 

The listed above values for the global characteristics of pp - interactions at the 
LHC are different from the other model predictions. First, the total cross-section 
is predicted to be twice as much the common predictions in the range 95-120 mb 
[|24ll25l due to strong increase of the a e i{s). The prediction for the inelastic cross- 
section is a ine i(s) = 75 mb and coincides with the predictions of the other models 
[|25l . Thus, reflective scattering would not result in worsening the background 
situation at the LHC, since the elastic scattering provides only two extra particles 
in the final state. However, the prediction for the total cross-section overshoots 
even the existing cosmic ray data. But extracting the total proton-proton cross 
sections from cosmic ray experiments is far from being straightforward (cf. e.g. 
Il24l 1251 and references therein). Indeed, those experiments are sensitive to the 
model dependent parameter called inelasticity and they do not the elastic channel. 
If we will use for the ratio a e i/a to t the value 0.67 at the LHC energy y/s = 14 
TeV and recalculate the total cross-section obtained from the cosmic ray data we 
will get the value about 227 mb in good agreement with the predicted value. This 
agreement, however, is merely an indirect confirmation of the model prediction 
for the total cross-section, since as it was noted there is an ambiguity in the elastic 
cross-section determination. The behavior of the ratio o e i/a to t at s — > oo does not 
imply decreasing energy dependence of a ine i. The inelastic cross-section <Ji ne i 
increases monotonically and grows as Ins at s — > oo. Such a dependence of 
a ine i is in good agreement with the experimental data and, in particular, with the 
observed falling slope of the depth of shower maximum distribution [26J. We will 
discuss some of the cosmic ray related issues in the next section. 

3 Reflective scattering in cosmic rays 

As it was already noted, the important role in the studies of cosmic rays belongs 
to the inelasticity parameter K, which is defined as ratio of the energy going to 
inelastic processes to the total energy. The energy dependence of K is not evident 
and cannot be directly measured. The number of models predict its decreasing 
energy dependence while other models insist on the increasing energy behavior at 
high energies [|27l . Adopting simple ansatz of geometric models where parameter 
of inelasticity is related to inelastic overlap function we can use the following 
equation [|28) 




to get a qualitative knowledge on the inelasticity energy dependence. 
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The estimation based on the particular model with antishadowing ll23l leads 
to increasing dependence with energy till E ~ (3 — 4) ■ 10 7 GeV. In this region 
inelasticity reaches maximum value (K) = 1, since a e i/a to t = 1/2 and then starts 
to decrease at the energies where this ratio goes beyond the black disk limit 1/2. 
Such qualitative non-monotonous energy dependence of inelasticity is the result 
of transition to the reflective scattering regime. 

Reflective scattering results in relative suppression of particle production at 
small impact parameters: 



due to the peripherality of da ine i/db 2 . Thus, the main contribution to the integral 
multiplicity n(s) comes from the region of b ~ R(s) and the distinctive feature 
of this mechanism is the ring-like shape of particle production which will lead 
to correlations in the transverse momentum of the secondary particles. It means 
that the enhancement of particle production at fixed impact distances b ~ R(s) 
would lead to higher probability of the circular events. Such events would reflect 
the production geometry with complete absorption at the impact distances equal 
to the effective interaction radius R(s). The enhancement of the peripheral par- 
ticle production would destroy the balance between orbital angular momentum 
in the initial and final states; most of the particles in the final state would carry 
out orbital angular momentum. To compensate this orbital momentum the spins 
of the secondary particles should become lined up, i.e. the spins of the particles 
in the ring-like events should demonstrate significant correlations. Of course the 
observation of such effects is difficult due to the multiple interaction of the sec- 
ondary particles in the atmosphere. However, the circular event has been observed 
experimentally ||29l . 

It might be useful to note that the rescattering processes are affected by the 
reflective scattering. At the energies and impact parameters where reflective scat- 
tering exists there will be no rescatterings, i.e. the head-on collisions will experi- 
ence a reflective elastic scattering at the energies s > sq and due to this fact less 
number of secondary particles will be detected in EAS at the ground level. This 
should provide a faster decrease of the energy spectrum reconstructed from EAS, 
i.e. it will result in the appearance of the knee. Thus, the hadron interaction and 
mechanism of particle generation will be changing in the region of ^fs = 3 — 6 
TeV. Indeed, the energy spectrum which follows a simple power-like law ~ E^ 
changes its slope in this energy region, i.e. index 7 increases from 2.7 to 3.1. The 
interpretation of the cosmic-ray data is complicated since the primary energies of 
cosmic particles are far beyond of the energies of modern accelerators with fixed 
targets and existing simulation programs merely extrapolate the present knowl- 
edge on the hadron interaction dynamics in the unexplored energy region. This 




(16) 
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might be an oversimplification and here we would like to interpret the existence 
of the knee in the cosmic ray energy spectrum as the effect of changing hadron in- 
teraction mechanism related to the appearance of the reflective scattering at small 
impact parameters. 

Conclusion 

In this note we considered saturation of the unitarity limit in head-on (and small 
impact parameter) hadron collisions, i.e. when S(s, &)|&=o — > — 1 at s — > oo. 
Approach to the full absorption in head-on collisions, in another words, the limit 
S(s, 6) 1 6=0 — ► at s — > oo does not follow from unitarity, it is merely a result of 
the assumed saturation of the black disk limit. This limit is a direct consequence 
of the exponential unitarization with an extra assumption on the pure imaginary 
nature of the phase shift. On the other hand, the reflective scattering is a natural in- 
terpretation of the unitarity saturation based on the optical concepts in high energy 
hadron scattering. Such reflective scattering can be interpreted as a result of the 
continuous increasing density of the scatterer with energy, i.e. when density goes 
beyond some critical value relevant for the black disk limit saturation, the scat- 
terer starts to acquire a reflective ability. Having in mind the quark-gluon hadron 
structure it would be tempting to find the particular microscopic mechanisms re- 
lated to the collective quark-gluon dynamics in head-on collisions which can be 
envisaged as the origin of the reflection phenomenon. One can try to speculate at 
this point and relate the appearance of the reflective scattering to the Color-Glass 
Condensate in QCD (cf. e.g. QUI and references therein) merely ascribing the 
reflective ability to Glazma. However, as it was mentioned in the Introduction, 
the evident obstacle which preclude a direct link to the microscopic mechanisms 
is the problem of confinement. The unitarity and consequently its saturation and 
reflective scattering are the concepts relevant for the hadronic degrees of freedom. 
The concept of reflective scattering itself is general, and results from the S-matrix 
unitarity saturation related to the necessity to provide an indefinite total cross sec- 
tion growth at s — > oo. 

Thus, at very high energies (s > s ) two separate regions of impact param- 
eter distances can be anticipated, namely the outer region of peripheral colli- 
sions where the scattering has a typical shadow origin, i.e. S(s, 6)|&>.r( s ) > 
and the inner region of central collisions where the scattering has a combined 
reflective and absorptive origin, S(s,b)\b<R( s ) < 0. The transition to the nega- 
tive values of S leads to the appearance of the real part of the phase shift, i.e. 

<fo(s, &)!&<«(*) = 7T/2. 

The generic geometric picture at fixed energy beyond the black disc limit can 
be described as a scattering off the partially reflective and partially absorptive disk 
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surrounded by the black ring which becomes grey at larger values of the impact 
parameter. The evolution with energy is characterized by increasing albedo due to 
the interrelated increase of reflection and decrease of absorption at small impact 
parameters. This picture implies that the scattering amplitude at the LHC energies 
is beyond the black disk limit at small impact parameters (elastic S'-matrix is 
negative) and it provides explanation for the regularities observed in cosmic rays 
studies. 
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